Temperature dependence of the gain profile for THz quantum cascade lasers 
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We study the rapid decrease of peak gain in resonant-phonon THz Quantum Cascade Lasers with 
increasing temperature. The effect of various microscopic scattering processes on the gain profile 
as a function of temperature is discussed. We argue that increased broadening, primarily due to 
increased impurity scattering, and not diminishing population inversion, is the main reason for the 
reduction of peak gain. 



The Quantum Cascade Laser (QCL) fT] has proven 
to be a compact and robust coherent light source in 
the mid-IR region. The double-phonon resonance de- 
sign has given the mid-IR QCL an operating range up 
to and above room temperature QCLs in the THz- 
region are still restricted to low operating temperatures 
which limits practical applications. To this date, opera- 
tion up to 164 K in pulsed and 117 K in cw mode have 
been demonstrated for THz QCLs [3]. The common ap- 
proach to stabilize population inversion by the double- 
phonon resonance has not been successful for THz laser 
[j]. Extending the operating region of these devices to- 
ward higher temperatures requires a quantitative under- 
standing of the various microscopic temperature effects. 
Previously, theoretical studies have focused on the pop- 
ulation inversion decrease with temperature [1| and par- 
asitic transport channels Q. Here we show that even 
for constant population inversion, the lasing performance 
can be strongly reduced with temperature in THz struc- 
tures due to increased broadening of the gain profile. 

In order to study these effects, the THz QCL 
from Ref. 0] has been simulated. It is a four- well 
GaAs/Alo.isGao.ssAs design, see Fig. [H^a), emitting at 
1.9 THz (7.9 meV). Lasing was demonstrated up to 95 K 
in continuous wave operation and to 110 K in pulsed 
mode. 

Our transport calculations are based on the Non- 
Equilibrium Green's Functions (NEGF) method [1, Q. 
This method can be seen as an extension of density ma- 
trix theory [l^, 11 1 to resolve broadening effects. We 
include electron-phonon, impurity, and interface rough- 
ness scattering [Ij 



13[ within the self-consistent Born 



approximation, and electron-electron interaction within 
the mean-field approximation. Alloy scattering is 
expected to be negligible for the GaAs/AlGaAs mate- 
rial system discussed here, since the well material is a 
binary compound. From the steady-state solution to the 
transport problem, see Fig. [Ijb), the gain is calculated 
according to the full theory in Ref. [l5[. In our model 
the temperature enters in two ways: (i) the occupation 
of phonon modes and (ii) in the Debye approximation 
for the screening of ionized dopants. Throughout this 
letter, we refer to the lattice temperature T, which is typ- 
ically larger than the heat-sink temperature There 
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FIG. 1: (a) Conduction band profile with the modulus square 
of the six lowest Wannier-Stark states per period at the design 
bias of 46.5 mV/period. For clarity, the lowest state (!') of the 
previous period is plotted as well. The lasing transition occurs 
between states 3 and 4. The barrier at z ~ 30 nm is 5-doped 
with a sheet density of n = 2.25 x 10^° cm~^/period. (b) 
Calculated current-voltage relation. The experimental data 
shows the same shape and trend However, the peak at 

35 mV/period at low temperatures, corresponding to resonant 
tunneling from the injector state to the lower laser state (1' - 
3), is not observed and also the experimental current is 50 % 
smaller, both at high and low bias. 



is a monotonic dependence between both temperatures, 
which can be established by a detailed thermal modeling 
17| , so that the general trends of performance are similar 
for both temperatures. 

The obtained gain spectra are depicted in Fig. [21 The 
gain peak at 9 meV is strongly reduced with increasing 
temperature, in agreement with the observed vanishing 
of lasing at 110 K. A similar behavior is found for the 
absorption peak at 14 meV, corresponding to absorption 
between states 4 and 5. 

Commonly, disappearing gain is related to vanishing 
population inversion. Fig. [3l[a) shows that the popula- 
tion difference between the lasing states (3 and 4) indeed 
decreases with temperature. The occupation in state 3, 
assuming thermal equilibrium with respect to the heav- 
ily occupied state 1, is also depicted. The similar shape 
of the two curves suggests that thermal backfilling essen- 
tially stands for the reduction in population inversion. 
Since the lower laser state is separated from state 1 by 
approximately the optical phonon energy, this effect be- 
comes relevant only at temperatures far above 100 K. 
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FIG. 2: Gain spectra for different temperatures and a bias 
drop of 46.5 mV/period. 
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FIG. 3: (a) Subband population of the two lasing states as 
a function of temperature. The population of the lower laser 
state assuming thermal equilibrium with state 1 is also plotted 
(dashed) (ni ~ 1.3 x 10^" cm~^/period). (b) An approximate 
linear reduction of population inversion with temperature is 
found, while the gain at the lasing energy drops more rapidly. 
Estimating the mirror losses, qm ~ 1-3 cm~\^ and the waveg- 
uide loss, aw ~ 4-5 cm~^, according to Ref. [31, the maximum 
lasing temperature is below 100 K at this operation bias of 
46.5 mV/period. 



Therefore, this decrease in population inversion is much 
less compared to the significant reduction of gain in the 
studied temperature range, see Fig. [3{b). 

The excess reduction in peak gain can be related to a 
change in line shape. Indeed, Fig. [5] shows that the gain 
transition is broadened with temperature, which reduces 
its peak value. The same holds for the absorption tran- 
sition, which starts overlapping with the gain region at 
higher temperatures, causing an additional reduction of 
gain at the lasing energy. 

The included temperature dependent scattering mech- 
anisms are as follows: 

(z) Optical phonon emission, the dominant scattering 
effect in most QCLs, is proportional to nB(£'opt) + 1, 
where ^^(iJopt) is the Boltzmann-factor at the optical 
phonon energy. The phonon absorption rate is propor- 
tional to UBiEopt)- At T = 100 K, ns(£;opt) ~ lO'^ for 



an optical phonon energy, i^opt, of 36.7 meV. Thus, spon- 
taneous optical phonon emission is by far the dominant 
phonon process, which is temperature independent. 

(ii) Acoustic phonon scattering contributes very lit- 
tle to the total scattering. Assuming elastic scattering 
and high temperature compared to the typical acoustic 
phonon, acoustic phonon scattering contributes less than 
1% to the total elastic scattering and can be neglected 
(see below). 

(m) Impurity scattering is strongly influenced by the 
screening of ionized dopants by electrons, which is a 
complex many-body problem where approximations are 
necessary 19*. A common simple, but robust, ap- 
proach is the Debye-approximation where the electrons 
contributing to screening are assumed to be in ther- 
mal equilibrium and obey Boltzmann statistics. This 
results in a temperature dependent screening length, 
'^Debye — ^s^okBTe/ e^uzD, whcrc TisD IS the average elec- 
tron concentration, and Cseo is the dielectric constant 
of the sample. The idea is that a hot electron gas is 
less affected by the impurity potential and will, there- 
fore, screen it less. The choice of the electron tempera- 
ture Te and the three-dimensional electron density n^o 
is an intricate issue. For typical parameters of the sam- 
ple considered here, the screening length is of the or- 
der of the period. Thus, electrons from all subbands 
are expected to contribute to the screening [2^. As 
experiments for a similar structure show that the ma- 
jority of electrons resides in subband 1, with an effec- 
tive temperature close to the lattice temperature [l6| . 
we approximate Tg = T and use the average electron 
density over the entire structure in usd. At low tem- 
peratures compared to the Fermi energy Ep screening is 
better described using the Thomas- Fermi approximation, 
A^p = 2eseoi?F/3e^n3D. In the calculations, the Thomas- 
Fermi model is used at temperatures below the crossing 
point of the screening lengths for the two models, i.e., 
ksT = 2EplZ and vice versa. 

We use a simple screening model in order to show the 
relevance of the mechanism. For a full quantitative de- 
scription of screening in the nonequilibrium multisub- 
band system of the QCL, a much more detailed study 
is desirable. Here, both aspects of the spatially inhomo- 
geneous electron distribution [19|] as well as the frequency 
dependency of the dielectric function [2l|, in particular 
close to the optical phonon and plasma frequency, are of 
relevance. In this context it is also questionable if the dis- 
tribution of the subbands are sufficiently well described 
by an effective temperature at all. All these aspect go, 
however, far beyond our present study. 

In order to quantify these scattering processes, the 7- 
factors, describing scattering matrix elements [9l.ll5|. are 
plotted in Fig. [D The 7-factors relate to the scattering 
rate from state a to /3 by Fq^^j = "faafsp/h. The sum of 
the diagonal 7-factors, representing intraband scattering, 
for different scattering processes is a good measure of the 
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FIG. 4: 7-factors for intra-subband impurity scattering for 
the upper (solid) and lower (dashed) laser states for the THz 
(a) and mid-IR device (b). The bright red curve in (a) shows 
the FWHM estimated for the gain peak in Fig.[2l In (b), the 
bright red curve displays 7'™^ for the injector state (1). The 
kink in the 7'™^ at 82 K is the cross-over from the Thomas- 
Fermi to the Debye screening model (10.7 K in the THz de- 
vice), which is smoothed out in a more detailed calculation. 



states are only weakly affected by impurity scattering. 
In contrast, states with strong spatial overlap with the 
dopants, e.g., the lowest injector state 1, exhibit large 
(and temperature dependent) broadening, see Fig.lDJb). 

In conclusion, our simulations show that lasing of THz 
QCLs can terminate at lattice temperatures where still a 
large portion of the population inversion remains. We at- 
tribute the excess gain reduction to increased scattering- 
induced broadening due to a reduction of screening with 
temperature. Thus, the common focus on population 
dynamics is likely to miss an additional mechanism of 
temperature dependent gain reduction. 

The authors thank S. Kumar and C. Weber for helpful 
discussions and gratefully acknowledge financial support 
from the Swedish Research Council (VR). 



energy width of a state. We observe a strong increase in 
the 7™P-factor with temperature for both lasing states. 

Assuming high temperature compared to the typical 
acoustic phonon energy and linear phonon dispersion, 
the 7-factor in the elastic scattering approximation for 
acoustical phonon scattering reads 



(1) 

where m* is the effective mass, D is the deformation po- 
tential, V is the sound velocity, p is the mass density, 
and ipaiz) is the wavefunction corresponding to Wannier- 
Stark state a. Using material parameters for GaAs from 
Ref. [iil, we obtain 71^33 = 0.024 meV at 100 K, sug- 
gesting that acoustic phonon scattering can indeed be 
neglected in this context. 

In the THz device, interface roughness scattering 
is small compared to impurity scattering, 7333!'' = 
0.12 mcV and 7444!^ — 0.15 meV, the same holds for 
the optical phonon emission j^^li — 0.45 meV from the 
lower laser state. Therefore, the FWHM shows a strong 
temperature dependence following the dominating impu- 
rity scattering, see Fig. (H^a). Note that the FWHM is 
smaller than the sum of the individual widths due to 
correlations in the scattering environment [lB | . 

In mid-IR devices however, due to large conduction 
band offset between well and barrier material, interface 
roughness scattering can be stronger than impurity scat- 
tering. For the device presented in Ref. 23 1, for which our 
results are shown in Fig. IHJb), we find 7gggg = 1.36 meV, 
7ggg|^ — 1.15 meV, and 7gggf^ = 11.7 meV, which clearly 
dominates the total widths of the upper and lower laser 
states. Therefore, no strong temperature dependence of 
the gain profile is expected. An important aspect here is 
that the doping is located in the injector and is strongly 
screened by the large electron density, so that the lasing 
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